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Abstract  -/■<  *  '  -  ^  /-'«**- 

'^This  paper  is  the  second  in  a  series  of  three  that  analyze  a  method  of 
dimensional  reduction.  It  contains  some  results  for  approximation  of  functions 

4) 

on  the  interval  [-1,1]  with  elements  from  the  nullspace  of  P  ,  N  ^  1  ,  where 
P  is  a  second  order  ordinary  differential  operator.  A  special  case  of  this  is 
approximation  by  polynomials. 

The  one-dimensional  results  are  used  as  a  tool  to  prove  similar  versions  in 
several  dimensions.  These  multi-dimensional  results  are  directly  related  to  the 
approximate  method  of  dimensional  reduction  that  was  introduced  in  [13],  and  they 
lead  to  statements  about  the  convergence  properties  of  this  approach. 

The  third  paper,  which  analyzes  the  adaptive  aspects  of  the  method,  is 


forthcoming. 


1. 


Introduct ion 


In  a  recent  paper,  [13],  we  introduced  the  concept  of  dimensionally  reduced 
solutions  to  an  elliptic  boundary  value  problem.  These  are  obtained  by  projecting 
(in  the  energy)  the  true  solution  of  the  boundary  value  problem  in  the  n+l-dimen- 
sional  domain  w  x  [-h,h]  onto  spaces  of  the  form 

.  N 

VM  -  {  £  W  (x)4>.(y/h)|  w.  arbitrary)  , 

J-0  2  2  2 

00 

where  {<j>^)j=Q  is  a  given  set  of  functions  on  [-1,1],  (x  are  coordinates  on  a> 
and  y  ranges  over  [-h,h]).  For  some  basic  ideas  behind  this  concept,  see  the 
introduction  to  [13].  In  that  paper  the  focus  was  on  the  right  selection  of  the 
4>  j  *  s .  It  was  shown  there  that  for  a  very  wide  class  of  problems  the  4>j's  should 
be  selected  such  that 


spanUj}^”1  *  N(Pk)  , 

where  P  is  a  second  order  differential  operator  intrinsic  to  the  elliptic  boundary 
value  problem. 

The  estimates  of  the  error  given  in  [13]  were  asymptotic  in  h  ■*  0  .  The  present 
paper,  which  was  already  announced  there,  treats  convergence  as  N  -*■  ®  for  a  fixed 
value  of  h  .  For  convenience  the  fixed  value  of  h  is  set  equal  to  1. 

If  the  bilinear  form  associated  with  the  elliptic  boundary  value  problem 
satisfies  some  kind  of  "inf-sup"  condition,  then  it  is  well  known  that  the  rate 
of  convergence  is  the  same  as  the  rate  of  approximation  (cf.[l]). 

The  results  proven  here  are  hence  formulated  as  approximation-theoretic 
estimates,  and  as  such  have  Interest  regardless  of  the  concept  of  dimensionally 


reduced  solutions. 


2  1 

The  results  are  all  concerning  approximation  in  the  L  -  and  H  -norms, 

i.e.  ideally  suited  for  second  order  problems.  This  is  not  crucial  and  similar 

results  can  also  be  obtained,  e.g.,  for  the  norms  (/1||4^||2dy  +  J1!  |Aly,2u|  ^dy)1^2 

-1  dy  -1 

introduced  in  [13].  (A  here  denotes  a  strictly  positive  definite  (unbounded) 
operator  in  a  Hilbert  space  H  ,  and  u  is  a  function  with  values  in  H  . ) 

For  reasons  of  convenience  the  approximation  results  are  formulated  without 
any  boundary  conditions.  Various  types  of  fixed  boundary  conditions  can  immediate¬ 
ly  be  included  based  on  the  present  proofs. 

Estimates  of  the  error  introduced  by  dimensional  reduction,  as  N  goes  to  », 
do  exist  in  the  literature  (cf.  [5,7]).  The  problems  considered  in  those  two 
papers  come  from  structural  mechanics.  The  elliptic  operators  have  constant 
coefficients,  i.e.  the  <j>j's  are  polynomials.  The  results  are  not  nearly  as 
strong  as  the  ones  established  here.  In  [7]  the  estimates  are  based  on  the  degree 
of  regularity  in  C  -spaces;  this  is  not  very  well  suited  to  the  regularity  pro¬ 
perties  of  solutions  to  elliptic  boundary  value  problems,  and  therefore  gives 
crude  estimates.  The  estimates  in  [5]  are  based  on  bounding  the  remainder  in  the 
N-th  order  Taylor  expansion.  The  estimates  are  very  crude  and  do  not  give  any 
indication  of  the  rate  of  approximation. 

We  now  give  a  short  review  of  the  contents  of  this  paper.  In  section  2  it 
00 

is  shown  that  the  set  U  N(Pk)  (N  denotes  the  nullspace)  is  dense  in  H^"  for 

k=1  d  d 

any  second  order  operator  P  »  b  a  ,  where  both  a  and  b  are  bounded  from 
above  and  away  from  0  .  This  is  the  obvious  generalization  of  the  fact  that  the 

polynomials  are  dense  in  ,  and  it  also  justifies  the  claim  that  the  dimensionally 
reduced  solutions  introduced  in  [13]  will  get  arbitrarily  close  to  the  true  solu¬ 
tion.  In  section  3  the  rate  of  approximation  using  functions  in  N(P^),  N  >_  1  , 
is  linked  to  the  regularity  of  u  in  spaces  of  the  type  P(pm)  .  This  general 
result  though  is  not  always  optimal,  as  shown,  e.g.,  by  Theorem  4.1.  Section  4 


3 


is  devoted  to  giving  a  necessary  and  sufficient  condition  for  a  certain  rate  of 
approximation  by  polynomials  (i.e. ,  the  case  where  the  operator  P  is  a  constant- 
coefficient  operator).  In  section  5  this  is  carried  over  to  results  in  several 
dimensions  —  directly  relating  to  the  concept  of  dimensional  reduction.  The  exam¬ 
ple  treated  in  section  6  is  of  the  same  type  as  the  numerical  examples  in  [13]. 
Finally  the  appendix  contains  the  proofs  of  several  results  about  the  eigenvalues 
and  eigenfunctions  for  two-point  boundary  value  problems,  as  used  in  sections  2 
and  3. 

Note:  Unless  otherwise  stated,  all  constants  denoted  by  capital  letters 
are  generic. 
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2.  A  Density  Result 

Let  a  and  b  be  two  functions  in  L°°([-l,l])  such  that  3  constants 

a  ,  b  with 
o  o 


0  <  aQ  <_  a(y) 

°  <  bo  b(y) 

By  P  we  denote  the  differential  operator 


,  d  d 
dy  a  dy 

P  is  considered  as  a  mapping  L^([-l,l])  3  P(P)  -*•  L^([-l,l])  .  N(Pk)  denotes 

the  null-space  of  the  operator  P  for  any  integer  k  1  .  It  is  easily  seen 
k  1 

that  W(P  )  C  H  ([-1,1])  .  The  first  theorem  in  this  section  proves  the  density 
of  a  certain  class  of  functions  associated  with  the  operator  P  . 


Theorem  2.1 


U  N(Pk)  is  dense  in  H^([-l,l]) 
k*l 


Proof 

-  Y[  1 

By  a  change  of  variables,  y'  =  I  ds 

-1 


operator  P  transforms  into 


and  multiplication  by  -1  the 


d  a_  d_ 
dy '  b  dy ' 


We  can  therefore  for  the  proof  of  this  theorem  assume  that  P 

-  a(y)  t-  where  a  satisfies:  3  a  constant  a  with  0  < 
dy  dy  o 

Define  the  operator  Q  by 

P(Q)  -  P(P)nH1([-l,l]>  and  Q-P  on  P(Q)  .  Let  fQ 


is  given  by 
a0  la(y)  . 


denote  the  function 


5 


£0W 


TO  ds  €  N<p> 

-1 


and  define  the  sequence 


,  ,  00 

fi  i*0 


by 


f±  =  Q_i  f0  €  N(Pi+1)  . 


0  <  Xo  £  1  .  .  f_  X^  <_  ^  |  •  •  denotes  the  eigenvalues  of  Q  (repeated 

OO 

according  to  multiplicity).  Let  ^um^m=o  an  orthonormal  basis  of  eigenfunctions, 
u^  corresponding  to  X^.  f^  can  then  be  expanded  as 


f 


0 


00 


I 


a  u 
a  m 


» 


and  with  this  notation 


f 


i 


m=*0 


a  X-i 
m  m 


u 

m 


We  now  proceed  to  prove  that  any  eigenfunction  u^  can  be  approximated  from 

00  k 

within  (j  N(P  )  .  The  proof  is  by  induction  in  m  ,  and  we  start  with  m  *  0  . 
k-1 

For  any  i  >  1  we  have  that 


I  luo  *  “o1  xo  fil  l\  1 

H 

cii^^-o1  ‘ivn  *2 


»o  iml  (y°o)2  (v^f 1  i  c  (v°o)  (V‘1)21 1  «j 
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where  we  have  used  Lemma  A. 3  to  guarantee  that  a.  ^  0  .  Since  £  af  <  00 

J-l  j 

and,  by  Lemma  A.l,  A^/A^  <  1  this  shows  that 


-1  ,1  -  _ 
a0  X0fi  "  u0 

u  €  (j  JV(Pk) 
k-1 


as  i  -*■  »  ,  in  H^C-l,!])  , 


-  denotes  the  closure  in  H  )  . 

Now  assume  it  has  been  proven  for  some  m  >  1  that 


{u  J™-1  C  U  W(pk) 
1  j  j-0  -  k-1  (  ] 


From  Lemma  A. 3  we  know  that  a  ^  0  ,  and  hence  for  any  i  >  1 

m  — 


00  . 

u-a^A^f,  -x.-T  /  a  /«  Wa  A  V  u. 

m  m  m  i  m,i  ^  j  m  jj  j 


where  x 


:  «  -  £  (a  /a,  \(  A  /A  A  u  €  U  W(p* 

m,i  j-0  V  J  m)\  m  3/  j  k=1 


due  to  the  induction  hypothesis.  As  before  the  H  -  norm  of  the  sum 


can  be  estimated  by 


(^■X.'wr  (I,  -iJ 


Because  of  the  facts  that 


(,L  -ft 


j-m+1 


<  «  and,  by  Lemma  A.l,  A  /A  <  1 

m  m+i 


this  shows  that 


•SJW-*  - 


•  -Me 

x.+aXf.  -*■  u  as 
m,  l  m  m  i  m 


i  -V  oo  , 


in  Hx( [-1,1])  ,  i.e. 


{u.}”c  U  W(Pk) 
J  J  k«l 


This  finishes  the  induction  proof,  and  we  conclude  that 


{u  }"  C  u  W(Pk) 
J  -1  k=l 


From  the  definition  of  Q  it  immediately  follows  that 


0(Q1/2)  =  H1(f-l,l]) 


f  1  °°  1/0 

and  since  iu^ Jj=q  is  complete  in  V(Q  )  this  proves  that 


hx( [—1,1])  E  u  w(pk) 

k-1 

Now  if  u€H^([-l,l])  we  shall,  by  choosing  c  -  u(-l)  and  d  = 
1  , 

(u(l)-u(-l ))//  a  x(s)ds  ,  obtain  that 
-1 


u  -  C  -  df0  €  H( [-1,1]) 

Since  1,  f^  €  N(P)  we  see  by  a  combination  of  this  and  the  previously 
proven  inclusion  that 


^([-1,1])  -  U  N(pk) 
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Based  on  Theorem  2.1  we  can  easily  prove  a  result  concerning  the  dimensionally 
reduced  solutions  as  introduced  in  [13].  This  result  guarantees  the  fulfilment 
of  the  goal  stating  that  the  dimensionally  reduced  solutions  shall  be  able  to  get 
arbitrarily  close  to  the  true  solution. 

Let  u  denote  a  domain  in  IR11  with  a  Lipschitz  boundary. 

Theorem  2 . 2 

The  set 


w  *1  1 

{  £  v.  (x)q.  (y)  |  J  €  IN,v,  6H  (oj)  and  q .  €.  ^  N(P  )  for  0  £  j  <_  J  } 

j-0  2  2  2  2  k-1 

is  dense  in  H^(w  x  [-1,1])  . 

Proof 

Follows  immediately  from  Theorem  2.1  and  the  fact  that 

J  1  1 
t  ^  v.  (x)w,  (y)  1 J  €  lN,v,  €  H  (oi)  and  w.€H  ([-1,1])  for  0  <_  j  £J} 

j«0  -1  2  2  2 

is  dense  in  H^(w  x  [-1,1])  . 

a 
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3.  Estimates  of  the  Rate  of  Approximation 

In  the  previous  section  we  proved  the  density  of  a  certain  class  of  functions 

In  this  section  we  shall  prove 
some  results  concerning  the  rate  of  approximation.  The  first  theorem  is  the 
following. 


associated  with  the  operator  P  =  b  ~  a  . 


Theorem  3.1 

2 

Assume  that  a,b€C  ([-1,1])  ,  and  let  m  be  an  integer  0  .  For  any 

e  >  0  there  exists  a  constant  C  such  that 

£ 


inf 

v€N(PN) 


V(Pm) 


Vn  >  1 


Note.  ||*|]  denotes  the  norm  | |Pm(* ) | |  0  +  | | * | |  „ 

-  V(Vm)  L2  L2 

One  can  of  course  combine  the  statement  in  Theorem  3.1  with  interpolation  by  the 

K-method  (cf.  [4]).  This  way  it  follows,  that  if  u  €  CL2 for  some 

s 

0  <  s  <  1  ,  then  for  any  e  >  0 


inf 


u-v 


N, 


0  <  C  N 
.  2  —  e 


-ms+e 


v  €  A/(P  ) 


(L2,0(Pm)) 


s,°° 


2 

The  smoothness  requirement  that  a,b€C  ([-1,1])  is  not  necessary;  as  it 

2 

immediately  will  follow  from  this  proof  we  only  need  that  a/b  is  a  C  -function. 
This  last  remark  applies  to  all  of  the  results  in  this  section. 


In  order  to  prove  Theorem  3.1  we  need  an  auxiliary  result  concerning  uniform 
approximation  by  polynomials.  This  result  can  be  found,  e.g.,  in  chapter  6  of  [6]. 


10 


Lemma  3.1 

Let  4>  be  a  function  in  C^([c,d]).  Define  $  by  $(t)  =  cost  +  , 


t  €  1 0 ,  it  ] . 


Let  r  be  a  non-negative  integer.  There  exists  a  constant  Cf  such  that  for 


any  <}>  with  $€Cr([0,ir])  the  following  estimate  holds 


inf|*-PN|Q<  Cr(N+l)“r|;|r  V»>0 

Px, 


The  infinum  here  is  taken  over  all  polynomials  PN  of  degree  £  N  .  | • | 

i  0  r 

and  |*j  denote  the  norms  in  C  ([c,d])  and  C  ([o,tt])  respectively. 

We  now  continue  with 
Proof  of  Theorem  3.1 

Like  in  the  proof  of  Theorem  2.1  we  may  also  here  assume  that  P  is  given  by 

P  =-  ^  a(y)  .  Let  ^  i  >  0  be  defined  as  in  that  same  proof.  For  any  set 

N 

of  coefficients  {c^}^_q  we  ^ave»  using  Lemma  A.  3,  that 


u  -  1  c±£i  =  l  (S ./a.  -  l  c.X~  )u  , 

i-0  1*0  3  3  3  i-0  3  3 


00 

where  £  B.u.  *-s  expansion  corresponding  to  u  .  If  by  p 

j=-0  3  3  N 

J  r  i 

polynomial  p  (x)  =  £  c.x  ,  then  the  above  can  be  rewritten  as 

N  i=0  1 


If  by  p„  we  denote  the 
N 


"  -  J0  cifi  '  WaJ  ■  P»(XJ1))UJ 


and  this  leads  to  the  equality 


iu  -  i  cifiii22“  i  ai(ei/ai 

i-0  1  L  j-0  3  3  3 


PN(Aj1))2 
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As  In  the  proof  of  Theorem  2.1  Q  denotes  the  restriction  of  P  to 

P(P)  (1  H^([-l,  1] )  .  Let  us  now  for  a  while  assume  that  u€P(Q^)  . 

—  2  00 

Choose  A  so  that  {A^  }^agC[0,A]  .  Define  a  sequence  of  functions 
$m€C°°([0,.A])  ,  1  <_  M,  with  the  following  properties 

“  Sj/0tj  0  -j  -  M_1 
*=0  on  [O.Ajj1] 

Let  $  denote  the  mapping 

$(t)  ■  y  (1-cost):  [0,it  ]  •+  fO, A]  , 

it  then  follows  from  Lemma  A.l  and  Lemma  A. 2  that 

|*‘1(A~^1)  -  <f"1(A"1)|  >  C/j2  for  any  j  >  1  . 

This  estimate  tells  us  that  it  is  possible  to  construct  the  <j>M's  such  that 

VM,  j  >_  1 


iV$(t))lr  1  Cr  SUP  Kj+1)2r»^ai 
0<j<M-l  J  J 


Now  since  u€0(Qk)  we  know  that  1 6j  |  <  C^j+l)  2^|  |Qku|  |  2  , 

L 

this  with  Lemma  A. 4  we  get 


Sj/Ojl  ±Ck(j+l) 


-2k+l,  |rtk. 


Qu 


and  combining 


i.e.,  we  have  for  any  r  >_  k  and  M  ^  1 


I  V$(t))  ir  1  CrM2(r"k)+1||Qku||  2 

L 

Because  of  Letnma  3.1  we  can  now,  for  any  r  ^  k  ,  M  ^  1  ,  find 
of  degree  <_  N  such  that 


'VPnIo  1  Cr(N+l)'rM2(r-k)+1i|Qku|j  2 

L 

We  now  go  back  to  estimate  the  right  hand  side  of  (1) 

co  M-l  OO 

I  ^(0  /«  -pJJaT1))2  <  l  a?((* m-pJ)(a:1))2  +  2  l  e2  + 

j=0  J  J  J  w  J  j=0  2  w  2  j=M  2 

OO 

+  2  l  a‘«*  m-pS)(X:1))2 

j-M  J  2 

The  first  and  the  third  sum  can  be  estimated  by 

Cr(N+l)~2rM4(r"k)+2  ^  a2^ | | Qku | |22  < 

<  C  (N+l)"2rM4(r"k)+2| |Qku| |2  , 

L 

the  second  by 


C 


k 


r  (j+l)"4k!  lQku|  i22  1  C.M‘4k+1 
j-M  L 


polynomial 
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In  summary  we  have  therefore  proven 


N 


ll»  “  I  c  f  ||  <  C  (<N+l)“^rM' 

i-0  L 


-2r„4(r-k)+2  +  M~4k+1^ j jpk^  |2 

L2 


for  any  r  ^  k  ,  N  ^  0  and  M  ^  1  .  Taking  M  *  [v^J]  +  1  ([•  ]  denotes  the 
integer  part)  this  estimate  gives 


(2) 


l«-  l  c  f  ||2  <  C  (N+l)'2k+1|jQku||2  , 

1=0  1  T  K  * z 


L 


k  2 

all  provided  that  u€V(Q  )  .  Let  11^  denote  the  L  -projection  onto  linear 
combinations  of  the  functions  f^ . . . ,  fN  . 

(2)  expresses  that 


Mnu|  |  2  1  Ck(N+l)‘k+1/2l  |qku|  I  2  , 

L  Li 


at  the  same  time  it  is  clear  that 


I  lu"nNul  I  2  -  I  luN  2 

L  L 

Applying  interpolation  by  the  K-method  we  get  for  any  0  £  m  k 


-(k-l/2)f 

I h-nNu| I  2  <  Ck(N+l)  ( |Q  u| I  2 

L  L 

Now  let  m  be  fixed  and  k  -►  ,  from  the  previous  inequality  we  then  get 


Ve  >  0  3  C 


such  that 


lu-l^ull  2  1  Ce(N+l)“nri'e|  lQmu|  I  2 
L  L 


provided  u£P(Qm)  .  If  we  only  know  that  u€P(Pm)  ,  then  choose  {g^  }™=^C 


N (P  )  such  that 


*s'\  -  p3'x“  £°r  y  ~  ± 1 


p8j  -o 


for  y  =  +  1 


and  any  i  4  j-1  ,  (this  is  obviously  possible) .  This  way 


u-  [  gi  €V(Qm)  and 

j-1  2 


l  gl€N(Pm) 
j-1  2 


From  (3)  and  (4)  it  now  follows  that 


I |u  -  n  (u  -  l  g  )  -  f  8j II  2  — 
w  j=l  J  j-1  1  L 


Ce(N+l) 


•  llQ^-  l  g  >||  2-Ce(N+irm+e!!Pmu! 
J-1  2  L 


2 


Since  the  image  under  11^  is  contained  in  N(P  )  ,  this  estimate  yields  the 
desired  result  for  N  ^  m  .  There  are  only  a  finite  number  of  N's  <  m  ,  and 
hence  the  result  can  be  obtained  for  all  N  by  possibly  increasing  CE  . 
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Based  on  Theorem  3.1  we  can  prove  the  following  result  concerning  approximation  in 
.  1 

the  H  -norm. 


Theorem  3.2 

2 

Assume  that  a,b€c  ([-1,1])  ,  and  let  m  be  an  integer  >_  1  .  For  any 


e  >  0  there  exists  a  constant  C  such  that 

£ 


inf  |  j  u-v  |  |  ,  C  N  ^+£  |  |  u  |  |  V  N  >_  1 

v  €  N  (PN)  H  '  *<P"> 


Proof 

From  Theorem  3.1  it  follows  that  there  exist  v^€N(?  )  such  that 


||Pu4bI|  2  <ce  N-"-1+e  Hull  . 

N  s  V(Pm ) 

Now  choose  v  €  A/(PN+^)  with  Pv„  ■  v„  ,  and  such  that  v„  =  u  for  y  *  +  1 

N  N  N  N  — 

It  then  follows  that 


!'*-vnH  j  <  cjlPu-Pv  1  i  2-c|iPu-v  II  2  <  C  jj~nrt'1+e  ||u|| 

N  H1  «  l*  N  E  P(Pm) 


We  can  also  easily  prove  a  result  relating  to  the  dimensionally  reduced  solu¬ 
tions  as  introduced  in  [13],  Let  <aCIRn  be  a  domain  with  a  Lipschitz  boundary, 
x  *  (x^,...,xn>  denotes  coordinates  in  w  and  y  ranges  over  [-1,1].  P  denotes 
b(y)  -jp  a(y)  jp  considered  as  an  operator  L2(wx  f-l.ll  )DP(fr)  ->■  L2(u  x  [-1,1])  . 


m 
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Theorem  3.3 

2 

Assume  that  a,b€C  ([-1,1])  and  let  m  be  an  integer  >_  1  .  Let  u  be 

an  element  of  L  (w  x  [-1,1])  with  . . u  ,  u€P(Pm)  .  Then  for  any 

X1  n 

e  >  0  there  exist  C£  (independent  of  u  )  such  that 


inf  ||u-vj| 

HX(m  x 


r_i 


<  c  l 


Here  V, 


N 

denotes  the  set  {£  w  (x)4>  (y)  |v,  CH1^)},  where  {$ ,  is  a  basis 

j=0  J  J  J  J  J=u 


for  N(Pk)  C.  H1 ( [—1,1]) 


Proof 


Let  vN  denote  the 


N 

orthogonal  projection  of  Pu  onto  (  [  w  (x)<t>.  (y)  jw  €L  (<*•) 

1=0  J  J  J 

-z  3v 

in  the  L  (ui  x  [-1,1])  inner  product.  Then  it  is  clear  that  t— ! -  is  the  IT 

1 

projection  of  f*  u  onto  the  same  subspace.  From  Theorem  3.1  we  immediately  get 


n 


l  I  I* ‘dr<u_vH)  I  I  2  +  ll^-^ll  2  1 

i-1  i  LZ(o)  x  [-1,1])  N  L2(o)  x  [-1,1]) 


<  C£  (N+l)  -nr*'1+e 


£>(?“) 


for  any  function  v*  €VN+3  with  Pv*  -  vN  (if  N  is  odd  such  a  v*  will  be 

contained  in  Vj^  *  but  this  is  not  necessarily  so  for  N  even) .  Now  choosing 
vj$  so  that  also  v*  ■  u  for  y  *  +  1  (this  is  obviously  possible),  it  follows 


that 


mmmM 


4.--  "V. 
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4.  The  Constant  Coefficient  Case  -  1  Dimensional  Results 

The  following  two  sections  are  devoted  to  the  case  where  the  operator  P 

has  constant  coefficients.  In  the  previous  section  we  estimated  the  rate  of 

approximation  for  general  P's  ,  but  the  estimates  established  there  do  not  have 

exact  inverse  counterparts  nor  are  they  always  optimal.  As  will  be  shown  in  this 

section  and  the  next,  the  question  of  approximation  rate  can  be  much  further 

clarified  when  P  is  a  constant  coefficient  operator.  We  start  with  an  analysis 

1c 

of  the  1  dimensional  problem.  The  space  N(?  )  ,  k  1  ,  consists  simply  of  all 

polynomials  of  degree  2k-l  .  Theorem  3.1  combined  with  interpolation  says 

that  if  u€Ht([-l,l])  ,  then  there  exist  polynomials  p^  ,  of  degree  N  ,  such 

that  I  I u— p„  1  1  _  <  C  (N+l)  t^+€  .  Under  the  present  simplified  circumstances 
N 1  '  i.  —  e 
L 

we  can  prove  a  better  result.  In  the  formulation  of  this  result  we  use  the  Besov 

spaces  B~  ,  t  >  0  ,  (cf.  [4])  ,  instead  of  the  ordinary  Sobolev  spaces  Ht  . 

For  an  interpretation  in  terms  of  the  spaces  H1"  use  the  inclusions  HC  C  b!:  C_ 

i. 

e  valid  for  any  t  >  0  ,  e  >  0  . 

Theorem  4.1 

Let  t  be  a  given  positive  number.  There  exists  a  constant  such  that 

for  any  u€B2  ,,,([-1,1])  one  can  find  a  sequence  of  polynomials 
degree  of  <  N  ,  with 

1 lu_pNl I  2  -  ct(N+1)~tl !ul I  t 

L  B2,® 

Note.  A  similar  theorem  is  also  valid  for  the  H^-norm.  The  estimate  here 
becomes  (for  t  >  1) 


|u-p 


N1 


,1± 


ct(w-i) 


-t+1 


lu 


2  ,c 


i  a 


The  rate  of  approximation  established  in  Theorem  4.1  is  optimal  in  the 
following  sense. 

Theorem  4,2 
2 

If  u€L  ([-1,1])  and  there  exist  a  constant  C  and  a  sequence  of  poly- 
00 

nomials  ^e8ree  PN  <_  N  ,  such  that 

j | u-PN | |  2  <_  C(N+1)  C  for  some  t  >  0  , 

L 

then  «€«5,JlocnB^  . 

Note.  Theorem  4.2  is  not  an  exact  inverse  of  Theorem  4.1  since  it  only 
t/2 

guarantees  that  u€B_  ([-1,1])  .  But  based  on  Theorem  4.2  we  conclude  that 

■Z,® 

for  a  general  type  function  in  b!|  ([-1,1])  we  cannot  expect  more  than  an  approxi 

mation  rate  of  (N+l)  t  . 

Theorem  4.2  is  optimal  in  that  one  can  find  u  such  that  | | u— pN | |  <_  C£(N+1) 

and  *  »  u^B2+®^loc  for  any  e  >  0  ^cf-  • 

The  proof  of  Theorem  4.1  is  very  simple,  based  on  transforming  u  into  a 
periodic  function  and  estimating  the  remainder  of  the  k'th  order  Fourier  expansion 
Details  can  be  found,  e.g.,  in  [2]. 

The  proof  of  Theorem  4.2  is  not  quite  as  simple.  The  cornerstone  is  the 
so-called  Bernsteins  inequality 


^dy^N  I  ®l2 


valid  for  any  polynomial  of  degree  <_  N  .  For  more  details  see  [2]  or  [10]. 

As  already  noted  Theorem  4.2,  although  optimal,  is  not  an  exact  inverse  of 
Theorem  4.1.  This  can  be  taken  as  evidence  that  the  standard  Sobolev  or  Besov 


spaces  are  not  very  good  for  expressing  the  kind  of  regularity  needed  for  a 


certain  rate  of  approximation  by  polynomials.  They  do  not  take  into  account  the 
well  known  fact,  already  noted  by  Timan  (cf.  [9]),  that  the  polynomials  have  a 


certain  ability  to  absorb  singularities  at  the  end  points  of  an  interval. 


d  2d 

Let  L  denote  the  operator  -  -jj—  (1-y  )  with  domain  of  definition 


P(L)  *  {u£L2([-1,1]>  [  Lu€L2([-1,1])} 


Now  introduce  the  Besov  spaces  n  ,  t  >  0  ,  by 


H*  =  (V(Lp);V(Lq))e  y 

where  p,q  are  two  integers  with  0  £  p  <  t  <  q  and  0  <  s  <  1  is  selected  so 
that  p(l-s)  +  qs  **  t  .  Because  of  Theorem  14.1  in  [8],  which  says  that 

„  -  P(L t  and  the  reiteration  theorem  on  p.50  of  [4] 
it  follows  that  modulo  equivalent  norms  HZ  is  independent  of  the  choice  of  p 
and  q  . 

We  are  now  in  a  position  to  characterize  completely  the  regularity  needed 
for  a  certain  order  of  approximation  by  polynomials. 

Theorem  4 . 3 

Let  t  be  a  positive  number.  For  any  u €  HZ  we  can  find  a  sequence  of 
00 

polynomials  ^Pjj^n»0  ’  t*le  degree  of  pN  <_  N  ,  such  that 

I  I U-PN I  I  2  -  lul  i  t 

L  H 

2 

On  the  other  hand  if  u€L  ([-1,1])  and  there  exists  a  constant  and  a 

oo 

sequence  of  polynomials  {Pfj^N*o  *  ^e8ree  Pfj  —  ^  *  suc^  that 
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I lu“PNi !  2  1  cu<N+1>“2t 

L 

then  u€Ht  and 


I lul I  t  1  C(CU+| |uj  j  2) 

H  L 

for  some  constant  C  independent  of  u  . 

Note.  Cy  here  is  not  generic,  it  is  the  same  constant  in  the  two  inequalities. 
Proof 

We  start  by  proving  the  direct  part.  It  is  well  known  that  the  eigenfunctions 
of  L  are  the  Legendre  polynomials  {£.  •  Also 


L(^)  -  kCk+l)^ 


Let  now  u  be  an  element  of  V(LP)  and  let  £  a  £  be  the  Legendre  series 

m~0  m  m 


for  u  .  Since  uCV(LP)  we  know  that  7  a^m^  <  ||u||2 

_ r,  m  n 


N 

l  v  -l  ,  then 
m  m 

tn-0 


U3”0 


P(LP) 


Define  p 


N 


|u-PH| 


m*N+l 


a2  < 
m  — 


(N+l) 


-4p 


m-N+1 


2  4p 
am 
m 


<  (N+l) 


-4pi 


I2 

V(Lp) 


i.e.,  I  I u-p „ |  |  -  <  (I«-l)"2p|  |uj  | 

N  L2  P(Lp) 

Interpolation  applied  to  this  gives  the  desired  result. 

We  now  turn  to  the  inverse.  Assume  that  there  exists  polynomials  pN  of 
degree  £  N  such  that 


r 


22 


lu-PNl!  2  ±VN+1> 

L 


-2t 


Define  A.  ■  P,  ,  A.  ■  P  -P  ,  m  >  1  . 
o  I’m  2m  2m-^-  — 


Then 


I  Kl  I  2  -  cu  +  I  lul  I  2  and 

L  L 


I  KM  ,  <  I  |u-P  JL+I  lu-p  „  <  C/C2  m  >  1 


,-2tm 


“ML2-  "“W 


2m-l>  *2  -  °t  u 


2 

Since  | | Lp  | |  ,  <  Cn  | |p  |  |  for  any  polynomial  p  of  degree  £  n  ,  it  follows 
n  L  n  L  n 

from  above  that  for  any  non-negative  integer  q 


IU\ll  2  -  Cq,(Cu  +  INI  2?  and 
Xj  l 


(  ( La/t  |  |  .  <  G  -c.-c  .22(<J_t)m  m  >  1 
1  1  m1  '^2  —  q  t  u  — 


Now  define  v,  =  V  fi  .  We  then  get 
k  n  m 

m=0 


W‘WJo(|KM2+ll*-V 


iC,.t(Cu+  Hull  2  +  Cu  l  22(«-C)") 
T  m*l 


22(q_t)k(C.  +  Mu||  ,) 

4  f  L 


u  l2 


provided  q  >  t  .  At  the  same  time 


!u-vkllL2  =  1 lu“P2kl 'l2  -  Cu2 


By  defining  =  2  q  we  therefore 


(||u-vk||  +St||vk||  -)<C  (cu+  Hull  )  , 

L  o  / 1  g  \  ’  L 


P(Lq) 


and  since  s,  -*■  0  for  k  «  this  proves  that 


u€(L  2;V(Lq))t  =  Hl 


Nl  t  i  C  (C  +  Mu||  2) 
Ht  q,t  u  ^ 


This  theorem  also  allows  a  version  formulated  by  using  the  spaces  V(L*~) 
instead  of  the  corresponding  Besov  spaces.  It  is  derived  from  the  inclusions 
0(Lt)  C_HC  C  0(Lt  e)  valid  for  any  t  >  0  ,  e  >  0  . 

A  theorem  similar  to  Theorem  4.3  but  concerning  approximation  in  the 
norm  can  be  derived  immediately  based  on  Theorem  4.3. 

For  practical  purposes,  in  determining  the  rate  of  approximation,  the 
following  characterization  of  P(Lq)  (cf.  [3j)  will  often  be  convenient: 


P(Lq)  -  {u€L2([-l,l])!u€Hq([-l,l])  ,  (l-y2)qu€H2q([-l,l]> } 
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Let  us  now  give  one  simple  example  that  shows  how  a  result  similar  to 
Theorem  4,3  can  be  established  also  in  a  case  with  non-constant  coefficient. 

Example  4.1 

Let  P  denote  the  operator  a  1  a  »  where  the  function  a  is 


given  by 


a(y) 


for  y  ^  0 
for  y  <  0 


with  a+  and  a^  being  two  positive  constants.  (This  is  the  operator  arising 
in  the  numerical  examples  of  [13].) 

It  is  not  difficult  to  see  that  the  following  set  of  functions  is  a  basis 
for  W(PN) : 


*o  "  1  *  *i  “  7GT  y 

!  f7 

^2k-l(t)dt  *  ^2k+l  a(y)  ^2k(t)dt 

-1  -1 


K2k 


1  <  k  <  N-l 


Here  denotes  the  Legendre  polynomial  or  order  k  .  Performing  the 

Gram-Schmidt  orthogonalization  on  the  set  <f>Q ,  •  •  •  »^2N-2  **2N-1  ^in  that 

sequence),  in  the  inner-product  <u,v>  =  u(y)v(y)a(y)dy  ,  we  end  up  with 

a  new  set  of  functions  • • • »^2N-2’  ^2N-1  *  ^k  ds  a  Piecew*se  polynomial 

of  degree  k  .  Let  L a  denote  the  operator  -a  "L(y)  a(y)(l-y  )  ,  it  is 

then  clear  that 


Lj- 


ay2k 


k 

I 

j-o 


'3,k  ^2j 


0  <  k  <  N-l 
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La^2k-1  *  ^  Xj,k  *2j-l  1  -  k  -  N 


Now  we  have  because  of  the  orthogonality  of  the  and  the  fact  that 


La  is  self-adjoint  in  L  ( [-1,1] ,a(y)dy) 


<LaWa  =  ^k’^Ya  =  °  for  J  <  k 


i.e. , 


Yk  =  Xk\  for  any  k 


It  immediately  follows  that  A^  =  k(k+l)  . 

oo  2  co 

It  also  follows,  since  is  dense  in  L  ([-1,1])  ,  that  {  A^, 'Yk-O 

is  a  complete  set  of  eigenvalues  and  eigenfunctions  for  La  . 

As  in  the  proof  of  Theorem  4.3  we  now  get  that 


inf 

v€N(PN) 


<  CN 


-2t 


u€(P(i.P);  P(l5)) 

51  a  3 1 


for  any  0  <_  p  <  t  <  q  ,  and  0  <  s  <  1  chosen  such  that  t  =  p(l-s)  +  qs  . 

In  summary,  we  have  found  a  singular  operator  L  that  characterizes  the  rate 

a 

N 

of  approximation  with  functions  in  W(P  )  the  same  way  that  the  Legendre 
operator  does  with  polynomials. 
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5.  The  Constant  Coefficient  Case  -  Dimensions  Higher  Than  1 

In  this  section  we  prove  a  result  relating  to  the  dimensionally  reduced 

solutions  introduced  in  [13].  We  give  a  characterization  of  the  regularity 

needed  for  a  certain  rate  of  approximation.  The  approximating  functions  are 
N  1 

of  the  form  £  wi(x)P^(y)  where  w.  €H  (tu)  and  p  is  a  polynomial  of 
J  J  J  J 

degree  j  ,  i.e.,  the  operator  P  has  constant  coefficients.  u  as  before 
denotes  a  domain  in  IRn  with  a  Lipschitz  boundary  and  y  ranges  over  [-1,1]. 

In  the  proof  of  the  main  result  in  this  section  the  following  lemma  will 
be  very  useful. 

Lemma  5.1 

Let  C_  Hq  be  two  Banach  spaces  with  norms  | [• [ and  | |* | Iq  respec¬ 
tively. 

Let  w=0  be  an  increasing  sequence  of  subspaces  of  and  let  3  be 

a  positive  number.  We  assume  that  the  following  implication  holds 

u€Hn  and  inf  |  | u-q  j  L  <  C  (N+l)”6  ,  Vn  >  0 
q€VN 

u€Hl  and  I  |u|  | x  <  C(Cu  +  |  |uj  |Q) 

for  some  C  independent  of  u.  (C^  here  is  not  generic,  it  is  the  same  con¬ 
stant  in  the  two  inequalities.) 

As  a  result  of  this  it  follows  that  for  any  0  <  0  <  1 

u€Hq  and  inf  ||u-q||0<C  (N+l)“96  ,  VN  >  0 
q€VN 

u€(H0;Hl)e,»  and  ltu|l9>00£C(Cu+  ||u||0) 
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for  some  C  independent  of  u  .  (As  before  is  not  generic.) 

Proof 

Let  0  <  0  <  1  and  assume  that  there  exists  a  sequence  of  elements 
qN^VN’  N  0  ,  such  that 


Define 


I |u-qN! |Q  1  Cu(N+l) 


-06 


*0  =  ql 


K  =  q  -q  .  m>l  , 

m  2m  2m_1  ~ 


then 


I"  -  I  *.llo  ■  ic„  2 

m»0  2 


-k0B 


At  the  same  time 


1*0 1  <0  -  Cu  +  I  lu' '0  and 


M  lo  i  I  lo  +  I  I  i  ce,ecu  2'm86  .  "  i  1 


That  is,  A.  €V  m  >  1  ,  and 

m  2m  ~ 


2 


From  the  first  implication  in  the  statement  of  this  theorem  it  follows  that 


|2m6<9-1\jl1icf)j6  <CU  +  ||u||0  ♦  ||2“6<6-1\||0) 


We  therefore  get 


I|4.l|iiceyc«+' ll"ll0’-2 


kB(l-e) 


m=0 


,-k6 


If  we  define  =  2  the  following  inequality  now  holds 


Sk9(||u-  Jo  *»ll°  +  •kUJo  *JI1)  iC8,B(Cu+  N«ll0> 


Since  s,  -*■  0  for  k  -*•  «  this  proves  that 


u€(H0;H1)Qj0o  with 


l"lle,-ice,6<cu+  ll»ll0> 


□ 

Let  us  introduce  the  spaces 

KR  -  {u€U(LR)|^u€P(LR)  i-l,...,n}  , 

d  2  d  o 

L  here  denotes  -  (1-y  )  considered  as  an  operator  L  (u>  x  [-1,1]) 3 

^  2 

V(L)  -*■  L  (w  x  [-1,1])  ,  and  R  is  a  non-negative  integer. 

N 

VN  denotes  the  space  {  £  w .  (x)p .  (y)  jw.  CH1^) }  where  p  is  a  poly- 

j«0  J  J  J  J 

nomial  of  degree  j  ,  j  0  .  We  are  then  able  to  give  the  following  character¬ 
ization  of  approximation  by  the  spaces  V„  in  the  H^-norm. 

N 


l  t 


Theorem  5.1 


1  R 

Let  a  be  a  given  positive  number.  If  u€(H  (w  x  [-1,1]);  *>a/R,»  for 
some  integer  R  >  a  and  R  >  2a/e  ,  where  e  is  a  positive  number,  then  there 
exists  a  constant  C  such  that 

inf  |  | u—  q {  |  <  C<m)"2a+£,  VN  >  0  . 

q€VN  ir(u>x  [-1,1]) 

On  the  other  hand  if  for  some  e  >  0  there  exist  a  constant  C  such  that 


then 


inf  |  |  u-q  |  |  1  C(N+1)  2“~G  ,  VN  >_  0  , 

qCVN  r(ux  [-1,1]) 


u€  n 

R€m  ,R>a 


(H1(o>  x  [-1,1])  ;KR) 


a/R,< 


Before  we  proceed  with  the  proof  of  Theorem  5.1,  let  us  state  a  corollary 
that  immediately  follows  from  this  theorem. 

Modulo  e  this  is  the  equivalent  of  Theorem  4.3  in  more  than  one  dimen¬ 
sion. 


Corollary  5.1 

Let  a  be  a  given  positive  integer.  If  u€  fl  (H1(o)  x  [-1,1])  ;KR)  . 

R€IN,R>ct  a/K’ 

then  for  any  e  >  0  there  exist  a  constant  C£  such  that 


inf  |  |  u-q  |  |  ,  <_C  (N+l)~2“+e  ,  N>  0  . 

q€VN  IT(u)  x  [-1,1]) 


On  the  other  hand  if  for  some  e  >  0  there  exists  a  constant  C  such  that 


-2a-e 


inf  l|u-q|lHl(a)  x  {_1M  <C(N+1)  ,  VN>  0 


ue  n  (h^xkh);^) 


R€IN  ,R>a 


Proof  of  Theorem  5 . 1 

Assume  that  u€K  and  R  is  an  integer  >  1  .  Let  v  denote  the 

N  N 

"  r\ 

orthogonal  projection  of  Lu  onto  {  £  w  (x)p  .  (y)  |w.  €  L^(oj)  }  in  the  Lz(cu  x 

H3  ]  3  ]  2 

[-1,1])  inner  product.  Then  it  is  clear  that  vN  is  the  L  projection 

^  3  1 

of  L  -5 —  u  onto  the  same  subspace.  From  Theorem  4.3  we  immediately  get 


n  ^ 


J,  |lL^T<“-VN)ll,2,  .  .  ...  +  llUu-VB)l'I2.  .  .  ... 

i=l  i  L  (a)  x  [-1,1])  L  (to  x  [-1,1]) 


<  CrN 


i-1  3xi  P(lR )  P(lR) 


for  any  function  v*€V^  with  Lv*  »  v^.  Now  choosing  v*  so  that  also  /  v*(x_,y)dy 

,1  _1 

BJ  u(x,y)dy  for  any  jc€w  (this  is  obviously  possible)  it  follows  that 


£  II 3x  (u-vapl  I  ?  +  H^nII  2(  r  1  m  +  ^3y(u“vN^  ‘.2.  .  . 

i*l  i  L  (a)  x  [-1,1])  L  (to  x  [-1,1])  7  L  (oj  x  [-1,. 


°r  n2_2r  u| 


+  I  |Uu-v*)  I  I  2  )  1 

1])  N  LZ(u>  x  [-1,1])/ 


«l  I  'W.  +  I  M  I  'W, 
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Using  interpolation  on  this  result  we  get  that 


2  a 

inf  |  |  u-q  |  |  <_  C  N-2ct+  R  | 

?€Vn  H  (oj  x  [-1,1]) 


where  lilullla/R  a>  denotes  the  norm  on  (H^(a>  x  [-l,l];KR)a^R  .  Since 

2a/e  <R  ,  i.e.  2a/R  <  e  ,  the  direct  part  of  this  theorem  immediately  follows. 

Let  us  now  give  a  proof  of  the  second  part  of  the  theorem.  If 

-2R-  Re 

inf  |[u-q  ]|  .  £  CN  ""®"  for  some  R  >  a  ,  then  as  in 

q£VN  HX(u)  x  [-1,1]) 

the  proof  of  Theorem  4.3  it  easily  follows  that 


u €V(lR)  and  u €V(LR)  , 

i.e.  u<=KR  .  If  we  apply  Lemma  5.1  with  »  KR  ,  -  H^(co  x  [-1,1])  and 

9  =  a/R  ,  we  then  get  that 


inf  1| u-q  ||  <  C(N+l)'2a_e 

q€VN  HX(o>  x  [-1,1]) 


implies 

u€(H1(o)  x  [-l,l]);KR)a/R>00  , 


for  any  integer  R  >  a  ,  i.e. 


u  €  D  (H1  (a>  x  [-1,1])  ;KR)  /r  „ 

R63»,R>a  '  * 


□ 


i  k 
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For  the  conclusion  of  this  section  let  us  give  a  simple  example  that  shows 
the  practical  usefulness  of  Theorem  5.1  (or  corollary  5.1). 

Example  5.1 

Let  m  be  the  interval  [0,1]  .  Let  y  be  a  positive  number  and  let 
(A., 9)  denote  polar  coordinates  around  the  point  (1,1).  We  then  consider  func¬ 
tions  of  the  type  u  ■  AY$  ( 9)  ,  where  <t>  is  an  element  of  C°°([0,ir/2])  . 

It  is  not  difficult  to  prove  that 

u€  0  (H1^  x  [-1,1])  ;KR)  , 

R€IN,R>t  T/R’" 

for  0  <  t  <_  y  ,  and  that  for  a  general  choice  of  4>  this  is  not  so  for  any 

T  >  y  (if  Y  ^IN ,  then  this  is  not  so  for  any  <f>  and  x  >  y  except  =»  0) .  By 

an  application  of  corollary  5.1  we  therefore  get  that 

inf  |  |u-q|  |  <  C  (N+l)-2Y+e,  VN  >  0  , 

q£VN  H  (u  x  [-1,1])  E 

for  any  e  >  0  ,  and  at  the  same  time  that  for  a  general  choice  of  <J>  (or  for 

any  <j>  ^  0  in  the  case  y^IN)  there  exist  no  t  >0  and  C£  such  that 

inf  I  | u-q  j  |  <  C  (N+l)"2Y_e,  VN  >  0  . 

q€VN  H(w  x  [-1,1]) 

y 

A  function  of  the  type  K  <f(9)  is  a  typical  example  of  a  corner-singularity 
as  arising  from  the  solution  of  an  elliptic  boundary  value  problem. 

Theorem  5.1  (or  corollary  5.1)  is  thus  well  suited  to  predict  the  optimal 
order  of  convergence  (modulo  e  )  that  one  can  in  general  expect  by  dimensional 
reduction  of  elliptic  boundary  value  problems. 
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A  result  like  this  could  not  have  been  obtained  by  using  the  a-priori 
knowledge  of  the  regularity  of  solutions  to  elliptic  boundary  value  problems 
in  terms  of  standard  Sobolev  spaces. 


i 
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6 .  A  Simple  Example  of  Dimensional  Reduction 
Let  us  consider  the  boundary  value  problem 


Au  *  0 

in 

]0,l[x]-l,l[ 

u  -  0 

for 

x  **  0,1 

3u  .. 

an  ■ 8(10 

for 

y  -  +  i 

(n  is  the  outward  normal) . 

From  C 133  it  follows  that  the  optimal  choice  of  basis  functions  for  dimen- 

o 

sional  reduction  in  this  case  is  the  polynomials.  denotes  the  set 

N  ,  ox 

{  2,  w  (x)  p.  (y)  |w.£  H  ([0,1])}  ,  where  p  is  a  polynomial  of  degree  j  . 

j«0  J  J  J  2 


Let  u^  denote  the  projection  of  u  onto  VN  in  the  inner  product 


BC*,*) 


1  1 
r 


0-1 


(l±  ii.iiiij 

\3y  3y  3x  3x  / 


dydx 


It  is  clear  that 


inf 

o 


q 


€  V 


N 


2 

H1([0,l]x[-l,l]) 


_<  C  BCu-i^.u-Ujj)  £ 


<  C 


inf 

o 


q€V. 


N 


2 

^([o.iW-i.i]) 


and  hence  that  the  energy  error  B(u-u^,  u-u^)  is  asymptotically  in  N  equivalent 

1  0 

to  the  square  of  the  distance  (in  H  )  from  u  to  VN 
If  g  has  the  Fourier  series 


g(x)  -  l  g,  sin  kir 
k-1  K 


x 
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Chen  ic  immediately  follows  chat  u  is  given  by 


u<x,y) 


00 


l 

k-1 


cosh(kiTy)  ^k 
sinh(kir)  ktr 


sin  kirx 


In  terms  of  regularity  of  u  it  is  not  difficult  to  prove  that  this  formula  leads 
to  the  following  three  results: 


i)  V<*>0:  l  g2  k2a  <  -  <b^u€H:3/Z+a([0>l]x[-l,l]) 
k-1  K 


,3/2+a . 


ii)  \/“  1  0:  l  g£  k2a  <  « 


k-1 


fuc  n  (HX( [0,l]x[-l,  1  ])  ;<R)fl/ 

.  RGN  ,R>0  6/R> 

I  with  0  =  a  +  1/2 


Hi)  Va>0,  e  >  0:  u€  fl  (^(tO.lJxI-l.lJ)  ;^R)fl/R  „ 

R€IN,R>0  0/R» 

with  0  -  a  +  1/2  +  e  . 


r  2  2a  < 

^  l  8V  k  < 


k=l 


We  consider  two  different  choices  for  g 


g(x)  «  ir/4 
g(x)  -  x(x-l) 


For  the  first  choice  of  g  it  follows  that 


r  2  20 

Z  gv  k  <  »  for  any  0  <  1/2  ,  and 
k-1  K 


m  oo 


similarly  for  the  second  choice 


2  20 

l  g.  k  <  «  for  any  0  <  5/2  ,  and 


k-1 


2  ,5 


I  4* 


k=l 


Corollary  5.1  together  with  the  regularity  results  ii)  and  iii)  now 
ensure  that 

In  the  case  g(x)  «  ir/4 

B(u-Ujj,  u-u^)  will  converge  to  zero  faster  than  N  ,  y  e  >  0  ,  but 

on  the  other  hand  slower  than  N  ^  e ,  >  0 

In  the  case  g(x)  =»  x(x-l) 

-12+e 

B(u-uN>  u-un)  will  converge  to  zero  faster  than  N  ,  ye  >  0  ,  but 

on  the  other  hand  slower  than  N  ^  Z  ,V  e  >  0 

Figures  1  and  2  show  the  actual  computed  values  of  B(u-u^,  u-u^)  as  a 
function  of  N  in  the  two  different  cases.  Note  that  the  asymptotic  rate  of 
convergence  is  obtained  already  for  a  fairly  small  number  of  polynomials. 

For  details  concerning  the  computation  of  the  u^'s  see  [13]. 

Instead  of  using  Corollary  5.1  to  obtain  information  about  the  rate  of 
convergence  we  could  have  used  the  regularity  result  i)  and  a  2-dimensional 
version  of  Theorem  4.1.  This  way  we  could  at  most  have  predicted  convergence 
of  the  order  of  N  ^+e  and  N  ^+c  respectively,  i.e.  only  half  the  actual 
convergence  rate. 

In  C 13 3  we  considered  the  same  boundary  value  problem  as  here,  only  it  was 
on  the  domain  [0,l]x[-h,h]  for  some  h  >  0  ,  and  not  on  [0,l]x[-l,l]  .  From 


the  computational  results  there  it  follows  that  for  a  fixed  N  ^  2  , 
h  h  2 

B(u-u„,  u-u  )  behaves  like  h  ,  h  0  ,  in  the  case  where  g(x)  *  ir/4 
"  N  N 

(  uj:  is  the  projection  of  u  onto  {  £  w.  (x)  P . (y/h)  |  w  €  Hx([o,l]) }  )  . 

N  j.Q  J  J  J 

Comparing  this  to  the  result  obtained  here  for  g(x)  =  tt/4  it  is  seen  that  using 

N  polynomials,  y€C-l,lJ,  is  in  some  sense  equivalent  to  having  a  domain  of 
2 

thickness  1/N  .  A  similar  feature  has  been  noticed  by  comparison  of  the 

standard  h-version  of  the  F.E.M.  with  the  so-called  p-version  (cf.  [2j) . 

In  this  example  we  used  slight  variations  of  the  approximation  results 
proved  in  sections  4  and  5,  namely  with  fixed  boundary  conditions  =0  at 
x  *  0,1  .  The  proofs  of  these  results  follow  immediately  from  the  proofs  of 
the  similar  results  with  no  boundary  conditions. 


In  sections  2  and  3,  we  used  some  results  concerning  the  eigenvalues  and 


eigenfunctions  of  the  boundary  value  problem 

d  d 

-  ■% —  a  -7—  u  -  Xu 
dy  dy 

u(-l)  =  u(l)  =  0 

a  here  is  a  function  in  L  ([-1,1])  such  that  3  a  constant  a  with 

-*  o 

0  <  aQ  <_  a(y)  .  From  the  theory  of  Sturm-Liouville  systems  it  immediately 
follows  that  the  eigenvalues  (repeated  according  to  multiplicity)  form  a  sequence: 


0  <  X  <  X.  < 

o  —  1  — 


<  X  <  X_, 
—  m  —  m+1 


with  +  ®  as  the  only  limit  point. 

Lemma  A.l 

With  notation  as  above 


Xm  t  for  ®  t  » 

m  m 

i.e.  the  eigenvalues  are  all  simple. 

Proof 

Assume  that  for  some  m  j*  m'  X^  "  •  This  means  that  the  eigenvalue 

X  *  Xffl(*  X^,)  has  multiplicity  >  1  .  Let  u  and  u  be  two  linearly  independent 

*v, 

eigenvectors  corresponding  to  X  ,  and  let  v  ■  cu  +  du  be  a  nontrivial  linear 

dv 

combination  with  the  property  that  a  ■7—  *  0  for  y  ■  -1  (such  one  obviously 


A- 2 


exists) .  The  function  v  is  then  a  solution  to  the  initial-value  problem 


d  d  . 

-  —  a  —  v  ■  Av 
dy  dy 


0 


in  [-1,1] 
for  y  -  -1 


and  because  of  the  uniqueness  of  solutions  to  this  problem,  it  follows  that 
v  *  0  .  Since  v  is  a  nontrivial  linear  combination  of  u  and  u  this  shows 
that  u  and  u  are  linearly  dependent.  We  therefore  have  arrived  at  a  contra¬ 
diction,  i.e.,  A  j4  A  ,  for  m  j4  m’ 
m  m 

□ 

It  is  well  known  that  there  exist  constants  0  <  and  0  <  such  that 
C. (m+1)2  <  A  <  C0 (m+1)2  . 

By  imposing  an  extra  smoothness  requirement  on  a  we  can  obtain  a  much  more 
detailed  statement. 


Lemma  A. 2 

If  a  £  C2([-l,l])  ,  then 

A  -  (tr/4)2  »(m+l)2  +  0(1) 

m  1 

f  _i  /? 

where  l  •  j  (a(y))  dy 
-1 


A  proof  of  Lemma  A. 2  is  found  in  chapter  4  of  [11],  and  shall  not  be 
repeated  here. 

oo 

Let  {u  }  n  denote  a  sequence  of  normalized  eigenfunctions,  u  corres- 
m  m*»u  m 

ponding  to  Aq  .  Let  f^  be  given  as  in  section  2,  namely 


Lemma  A. 3 


The  function  fg  has  the  expansion 


00 


where  a  ^0  for  every  m  . 
m 


Proof 

That  fg  has  a  unique  expansion  is  wellknown.  The  coefficient 
given  by 


1  y 

r  r  i 

“m  "  j  j  7G)d3  Um(y)  dy 

-1  -1 


Now  assume  that  for  some  value  of  m  *  m  a  *  0  ,  i.e. 

o  m 

o 


i  y. 


-i  -i 


— 7— .  ds  u  (y)dy  *  0 

a(s)  m0 


Since  u  (y) 
m 

o 


1  d  d 

-  —  —  a  3—  u  we  get  that 

A  dy  dy  m 

m  J  J  0 

o 


;;  y 

-1-1 


^I)d9  [37  *  I7  U»J  <y)  dy  ' 0 


Performing  an  integration  by  parts  this  yields 


A-4 


and  the  last  integral  here  vanishes  due  to  the  fact  that 
We  therefore  conclude  that 


u  (1) 

IQ 

O 


U  (-1) 
m 


0 


d 

a  -r~  u 
dy  m 


u 

m 

o 


0 


for  y  >*  1 


On  the  other  hand  u  satisfies  the  differential  equation 

o 


d_ 

dy 


d 

a  —  u 
dy  m 


+  A 

m 

o 


in  C-lflJ 


Because  of  uniqueness  of  solutions  to  the  initial-value  problem,  this  implies 


that  u  *  0  .  We  have  thus  arrived  at  a  contradiction,  meaning  that  a  i*  0 
in  m 

o 

for  every  m  . 


□ 


Again,  by  imposing  an  extra  smoothness  requirement  on  a  ,  we  obtain  a 
very  detailed  result  concerning  the  decay-properties  of  the  am's 

Lemma  A. 4 
2 

If  a€ C  ([-1,1]),  then  ^  constants  0  <  and  0  <  such  that 


for  all  m 


Proof 

From  [ll]  p.  176  we  get  the  following  asymptotic  formula  for 


um(y) 


um(y) 


-  sr  «)] 


0((nrt-l)~2) 


where  5(y)  -  j  (a(s))  ds  and  l  *  C(l)  “  j  (a (a))  ds  .  The  function 


-1 


-1 


T  is  in  C  and  the  constants  D  satisfy 

m 


2  D  (independent  of  m  )  such  that 


1/D  <  ID  <  D  for  all  m 
in  — 


Also,  O(')  here  means  uniformly  in  y  .  Let  us  now  calculate  a 


m 


1  X 


m 


-1 


r  x 

f0<y)  um(y)dy  -  J  ^ds  um(y)dy  = 


-1-1 


.-2s 


X1  +  !2  +  0( (m+1)  *) 


I2  denotes  the  integral 


m+T  °m  .  l  j  ^)dsj  (a(y))"1/4  T«)cos  5  )dy  . 


By  a  change  of  variables  from  y  to  5  and  an  integration  by  parts,  it 


.-2s 


immediately  follows  that  ^  is  0((m+l)  ) 


1^  is  given  as 


1  /  X 


m 


dsl  (a(y))‘1/4  sin 


-1  \-l 


(*p«) 


dy 


By  a  change  of  variables  from  y  to  £  and  an  integration  by  parts  we 


IP 


r* 


A- 6 


I 


1 


l  (a(l))1/4  • 


1 

(m+l)u 


*(-l)m+1  +  0((m+l)'2)  . 


This  immediately  implies  the  existence  of  two  constants  0  <  and  0  < 
such  that 


m+1 


for  m  sufficiently  large.  Now  combining  with  Lemma  A. 3  and  possibly  changing 
the  constants  and  C2  we  get  the  desired  result. 
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